RIEMANN-HILBERT METHODS IN THE THEORY OF 
ORTHOGONAL POLYNOMIALS 

PERCY DEIFT 

To Barry Simon, on his 60th birthday. 
Mathematician extraordinaire, teacher and friend. 



Abstract. In this paper we describe various applications of the Riemann- 
Hilbert method to the theory of orthogonal polynomials on the line and on the 
circle. 



1. Introduction 

In this paper d/i denotes either a Borel measure on M with finite moments 

\xr dfi{x) < oo , m = 0,1,2,... (1) 

or a finite Borel measure on the unit circle 

dti{0) < oo. (2) 

SI 

In addition, unless stated explicitly otherwise, we will always assume that dfj, is a 
nontrivial probability measure, i.e. supp((i/Lt) is infinite and the integral of d/Lt is 1. 
Let 

p,,{x) = fc„a;" + . . . , fc„ > , n = 0, 1, 2, . . . (3) 

(/)„(z) = K„Z" + . . . , K„ > , 71 = 0, 1, 2, . . . (4) 

denote the orthonormal polynomials (OP's) with respect to d/i on R and respec- 
tively (see |Szej), 



Pn{x)pm{x)dp{x) ^ / (l3n{e ) 4'mie'' )dfl{9) ^ 5n,m , n, m > 0. (5) 

The fact that dfi is nontrivial implies, in particular, that the p„'s, and the c^n's, 
exist and are unique for all n > 0. 

As is well known, the p„'s satisfy a three-term recurrence relation 

lPn-l{x) + {x) = xpn{x) , n > (6) 

where 

a„ e M , 6„ > , n > (7) 
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and b-i = 0. Similarly the 0n's satisfy the Szego recurrence relation 
a/1 - |q!„P 4'n+iiz) = z(/)„(z) - a„(/)* (z) , n>0 

where 

|a„| < 1 , 71 > 



(8) 
(9) 



and for any polynomial q(z) of degree n 



q*iz) = z^ q{l/z) (10) 

denotes the so-called reverse polynomial. Following |Sim2| . we call the a^s Verhlun- 
sky coefficients. A simple computation shows that 

1 



Kn+l 



-0„+i(O) , n>0 



On M wc define the (rt + 1) x (n + 1) Hankel determinant 



D„ — det 



71 > , 



0<j,fe<n 

and on we similarly define the (n + 1) x (n + 1) Toeplitz determinant 



An = dot 



dti{9) 



n>0. 



(11) 



(12) 



(13) 



/ Q<j,k<n 

The determinants £)„ and A„ are closely related to the OP's {Pn},{4'n} respec- 
tively: Indeed one has (see for example \Sze\ ) 



Dn-l 



n > 1 



(14) 



Given dfi, the study of the algebraic and asymptotic properties of the quantities 
and also 

D„ and A„ , 

constitutes the core of the classical theory of orthogonal polynomials. 
The three-term relation ^ can be re-written in the form 

Lp{z) = zp{z) , p{z) = {po{z),pi{z),p2{z), . . .)^ , 

where L is an infinite Jacobi matrix, i.e. L is symmetric and tridiagonal 

/ao bo \ 
bo ai bi 

bi 02 ' ■ 

\o ■./ 

with > 0, i > 0. In the case that dfi has compact support on M, the operator L 
is bounded on 

£^ = i^u= {uo,ui,...)'^ : < oo| ■ 



(15) 



(16) 



Let 



F : {dfi on M : supp(ci/x) compact} {bounded Jacobi matrices on ^^}. 
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denote the map taking dfj, i—^ L. Conversely, if i is a bounded Jacobi matrix, then 
in particular L is self-adjoint, and we let d/i denote the spectral measure associated 
with L in the cyclic subspace generated by L and cq, where eo = (1, 0, 0, . . .)'^ G £2 ■ 
Thus 

1 \ /■ dii{x) 



eo,^eoj=y^, AeC\M (17) 

and it follows further that d/i has compact support. Let 

F : {bounded Jacobi matrices on l^} ~* {d-l^ on R : supp((i/i) compact} 

denote the map taking L to dfj,. The basic fact of the matter (see, for example, 
\K\. |Siml| . and also 1021 ) is that F and F are inverse to each other, F o F = id, 
F o F = id. From this point of view the (classical) orthogonal polynomial problem 
is the inverse spectral component of a spectral/inverse spectral problem. If the 
support of d/i is not compact, then the situation is similar, but the relation between 
dfi and L is more complicated because L is now an unbounded operator and we 
must distinguish between different self-adjoint extensions of L (see rA:, ^SiinTl for 
more details). 

In the case of measures d/i on the unit circle, the role of the Jacobi matrices is 
played by so-called CMV matrices C (see |Sim2| '). Such matrices C are unitary in 
£2 and pentadiagonal, and have the form 

C = LM (18) 

where L and M are block diagonal 

L = diag(eo,e2,e4,...) , M = diag(l, 61,63,...) (19) 

with 

Here 
and 

P, = ^l-\a,\^. (22) 

CMV matrices are named for Cantero, Moral and Velazquez |CMV| . but in fact 
they appeared earlier in the literature (see, in particular, IWatl l. Let 

^ : {dfi on S^} {CMV matrices} 

denote the map taking d/i — > C, the CMV matrix constructed from the Vcrblunsky 
coefficients aj = aj{dfi), j > 0, of d^, according to ((TH|l . ((T^ and Conversely, 
given a CMV matrix C, let d/i be the spectral measure associated with C in the 
cyclic subspace generated by C, C* — C^^ and eo. Let 

ip : {CMV matrices} {d^i on S^} 

denote the map taking C to d^. Then, as above (see jSim2| '). ip and ip are inverse 
to each other, and we see again that the classical orthogonal polynomial problem 
on is the inverse spectral component of a spectral/inverse spectral problem. 

The techniques used to analyze the direct spectral maps, F and ip, are generally 
very different from the techniques used to analyze the inverse spectral maps, F or 
■0, though sometimes there is some overlap (see e.g. |DK| '). It is also interesting to 
note that in the solution of integrable systems one needs knowledge of both F and 



) - J>0. (20) 



< 1 , J > (21) 
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F (or i/) and i/')- For example, the Toda lattice induces a flow Lo 
Jacobi matrices f |F)) 



L{t = 0)=Lo 



L{t) on 
(23) 



where 



60 ai &i 



61 02 



V 



\ 



B{L) 



I 



V 







61 



\ 



and the solution of H23() is given by the following well-known procedure ( M ): 



„2\t 



Lq dfiQ = F{Lq) dfitiX) 



L{t) = F{dfit 



The analysis of F and tp has benefited greatly from the powerful developments that 
have taken place over many years in the spectral theory of Schrodinger operators 
and their discrete analogs, reaching, over the last 20 years or so, and in the case of 
one dimension, a state of great precision. Here Barry Simon and his school have 
played a decisive role, and we refer the reader to jSim2| . in particular. Part 2. The 
systematic analysis of F begins with the classic memoir of Stieltjes 1894-1895. Up 
till that point, a great deal of information had been obtained concerning particular 
polynomials, such as Legendre polynomials, Jacobi polynomials, Hermite polyno- 
mials, etc., but a unified point of view based on the orthogonality relation (j^l had 
not yet emerged. The analysis of began in 1920, when Szego initiated the sys- 
tematic study of polynomials orthogonal with respect to a measure on S^, as in ||SJ). 
Szego's work in turn has led to many remarkable developments by researchers from 
all over the world, particularly the former USSR, Europe and the USA. We refer 
the reader to Simon's book |Sim2| . where these developments are discussed in great 
detail together with many fascinating anecdotes concerning their discovery. Start- 
ing in the early 1950's with the celebrated work of Gel'fand and Levitan, various 
techniques were developed to recover one-dimensional Schrodinger operators from 
their spectral measures. In the 1970's, techniques based on the inverse-Schrodinger 
method (see |CGe| and |F]) started to play a role in the analysis of F and V'- The 
goal of this paper is to describe one of these techniques, which is different from 
the techniques in |CGe| or 0, and which has proved extremely fruitful, viz., the 
Riemann-Hilbert (RH) method, also referred to as the Riemann-Hilbcrt Problem 
(RHP). The scope of the paper is limited to describing results for F and ip ob- 
tained by RHP. Some of the results that we describe are quite standard and are 
included only for purposes of illustration. Other results, particularly asymptotic 
results, have been obtained, so far, only through RH methods. For a full up-to-date 
discussion of what is known about F and ■0, including the seminal contributions of 
Golinskii, Ismail, Khrushchev, Lubinsky, Nevai, Rakhmanov, Saff, Totik and many 
others, we again refer the reader to jSiml| and |Sim2| . 

To begin, let S be an oriented contour in the complex plane C (see Figure^. 
By convention, if we move along the contour in the direction of the orientation, 
the (-l-)-side (resp. (— )-side) of the contour lies to the left (resp. right) (see again 
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Figure 1 . The contour S 

Figure 1). A k x k jump matrix t; on S is a mapping from S — !■ G£{k, C) such that 
V, G L°°(S). We say that an £ x fc-valued matrix function m{z) is a solution of 
the RHP (E, v) if 

(a) rn(z) is analytic in C\S 

(b) rn_|_(z) — m_(z)v(z), z G S, where m±{z) — lim m,{z') 

z' — '2:,2:G(±)— side 

If in addition i — k and 

(c) m(z) ^ / as z ^ oo 

we say that m is a solution of the normalized RHP 

Many technical issues arise. For example, in what sense do the limits m± exist? 
In what sense does m{z) ^ I in (c)? How should one understand (b) at points of 
self-intersection in S? Under what assumptions on S and v does a solution m{z) 
exist? And if we normalize as in (c), is the solution unique? We will not consider 
such issues here and in the text that follows, and we simply refer the reader to |CG| 
and the references therein for a general discussion of RHP's (see also |DZ4| for 
more recent information, and |BDTj for a discussion of points of self- intersections). 
In this paper we will consider almost exclusively problems with solutions m{z) 
that are analytic in C\I] and continuous up to the boundary and at 2 = oo. For 
such solutions, the limits in (b) and (c) are taken pointwise. Furthermore, for the 
problems we consider, the solution of the normalized RHP will always exist and be 
unique. 

At the analytical level, a normalized RHP is equivalent to a problem for coupled 
singular integral equations on E. This is seen as follows. 
Let denote the Cauchy operator on S, 

C^hiz) ^ j ^ z e C\I] (24) 

Jj2 s — z ZTri 

with boundary values 

{C^h){z)= lim {C^h){z') , zeS. (25) 

side 

Under reasonable conditions on E, g £(LP(S), the bounded operators from 
LP(S) LP(S), for any 1 < p < 00, and we have the relation 

Cf-C^^l. (26) 
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Let 

v{z) = {v-{z))-^v+{z) , zeS (27) 
be any pointwise factorization of v where 

v±, {v±)-' e . (28) 

Set 

w+=v+-I w^=I-v. 

and define the singular integral operator on S 

C^h = C^ihw^) + C^{hw+) (30) 

for row fc-vectors h. As w± £ i°°(S), e £(LP(E)), 1 < p < oo. Suppose in 
addition that 

w± e L^(S) for some 1 < p < oo , (31) 
and consider the equation for a fc x /c-matrix function /i 

(1 ~ CS)m = / (32) 

in / + LP(E), or more precisely 

(1 - cS)i. = C^i = c>_ + c^w+ e (33) 

where 

li = I + v , vGLP. (34) 
If a solution /i = / + of (|32|l - (|34|l exists, set 

m{z) ^ I + C^{fi{w+ + w^)){z). (35) 

Then a simple calculation shows that m± = and hence m+ = m_w, and as 

to(2) ^ / as z ^ oo, we see that (I35|l gives a solution of the normalized RHP 
(S,w). Thus the normalized RHP (S,w) reduces to the analysis of the singular 
integral equations (j221>- 

The connection between the OP problem and the RHP is due to Fokas, Its and 
Kitaev |FIK) . Let 

Pn = -^Pn = a:" + . . . , n > (36) 

denote the monic orthogonal polynomials associated with a measure 

dii{x) — w{x)dx , w{x) > (37) 

absolutely continuous with respect to Lebesgue measure on M, with x^w{x) G 
H^{M.), the first Sobolev space, for all j > 0. Let E = R, oriented from — oo 
to +CX), and equipped with jump matrix 

V = v{x) = ^ Q ^^i^ ^ ■ < X < oo. (38) 

Finally, for any n > 0, let X^"'> = iX^f)i<;j<2 solve the RHP {R,v) 

X("'(z) analytic in C\M 

X|"^(z) = X^_^'\z)v{z) , zeR (39) 



normalized so that 



V ,n / as z^oo 
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Then f |FIKj . in addition see |p2j) direct computation shows that 

X(^)(A-( -'"(^^ C{PnW){z) \ , . 

where C — denotes the Cauchy operator on E = R. In particular, 

P^z) ^ X['l\z) . (41) 

Furthermore, if Xl' denotes the residue of X'") ( q n j ^-t infinity, 



then 

and in the notation of ® 



ei = -^(^l"^)2i (42) 



a„ = {X[%^ - (Xi"+^^)n (43) 



= {XnMXr'%i (44) 

Also by lO and (|^ . 

-d7 = "^^""^ ^''^ 
Thus all the basic quantities of interest in the OP problem can be read off from the 
solution of the RHP (M,^) above. 

On the unit circle, the situation is similar. Let 

$„ = — 0„ = z" + . . . , n > (46) 

denote the monic orthogonal polynomials associated with a measure 

d^{e)^u;{e)^ (47) 

absolutely continuous with respect to Lebesgue measure on with uj{9) E H^{S^), 
uj{9) = uj{9 + 2tt). Fix 71 > and let S = S*^, oriented counterclockwise. Equip 
with the jump matrix 

V = m . ( I -(^f " ) , z = e- (48) 



and let F^") = (r/"^)i<.j<2 solve the RHP {S^ 



V 



• y(")(z) analytic in C \ 5"^ (49) 
. yI''\z) = Y}_''\z)v{e) , z = e''^ e (50) 

normalized so that 

• r(")(z)(^ ) ^/ as z^oo. (51) 

Then again (cf. |BDJ| ) direct computation shows that 



8 



PERCY DEIFT 



where C = C^^ denotes the Cauchy operator on and is the reverse poly- 

nomial as in 1)10(1 . In particular 

$„(z) = rA")(z) (53) 

and hence by (fTT|l . 

a„_i = -r/i"^(z = 0) . (54) 

Also 

«Li = ->^ir^(^-0) (55) 

and hence 

^ = -4"^(z = 0). (56) 

Again we see that all basic quantities in the OP problem on the circle are expressed 
in terms of the solution F^") of the RHP (5\i;). 

The outline of the paper is as follows. In Section 2 we show how to use the 
RHP's (R, v) and {S^,v) above to derive various identities, equations and formulae 
for the OP problem. In Section 3 we describe the application of the steepest de- 
scent method of Deift-Zhou for RHP's to asymptotic problems for OP's. Finally, 
in Section 4 we describe the application of RH ideas to areas related to the OP 
problem, such as random matrix theory, multi-orthogonal polynomials, orthogonal 
Laurent polynomials, and the rarefaction problem for the Toda lattice. 
Technical Remark In most of the paper we will be considering probability mea- 
sures with some degree of smoothness as in (|37ll and (|47|l above. For such weights 
we then use the RHP's to derive, in particular, various identities such as 
(|85|l . etc. If dfi{x) is an arbitrary probability measure on M with finite moments, or 
d^lO) is a probability measure on 5^, we can approximate d^{x) and d^{0) appro- 
priately with smooth measures dfi^{x) and dfi^{9) respectively: For such measures 
©, ®i l|85|l . etc., are true, and letting e J, we conclude that these identities 
are true, as they should be, for all measures dfj,{x) and dfi{6) as above. Similar 
considerations apply at many points in the paper and we leave the details to the 
interested reader. 

2. Applications of (M, w) and {S^,v): identities, equations and 

formulae 

The applications of Riemann-Hilbert techniques to OP's are principally of two 
types: 

(a) algebraic 

(b) asymptotic. 

Under (a), the goal is to derive identities, equations and useful formulae for the 
OP problems. Under (b), the goal is to determine the asymptotic behavior of the 
OP's Pn, Pm 4>m as n — !■ cx): Here one considers the case where the weight uj{x) 
is independent of n, as well as the case where to^x) depends on n in a prescribed 
fashion (see (|106|l below). We consider (a) in this section, and (b) in the next. 

Regarding (a), there is a general methodology, which may be traced all the way 
back to the original work of Gel'fand and Levitan, and which may be stated loosely 
as follows: If the jump matrix for a RHP is independent of a parameter, then 
differentiation with respect to that parameter (or taking differences in the discrete 
case) leads to an equation/identity. 
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We illustrate this methodology, which may be viewed as the analog for RHP's 
of the celebrated theorem of Noether on conserved quantities for dynamical sys- 
tems, first in the case of the defocusing Nonlinear Schrodinger Equation (NLS). In 
1975 Shabat observed that the inverse scattering problem for the one-dimensional 
Schrodinger equation could be rephrased as a RHP. Because of the connection 
between Schrodinger operators and the Korteweg-de Vries (KdV) equation, this 
meant that KdV, and by extension all l-|-l-dimensional integrable systems, could 
be solved by a RHP. In the case of defocusing NLS, Shabat 's observation amounts 
to the following (see e.g. |DZ4p . Let q{x, t) be the solution of NLS on the line 

iqt + qxx-'2\q\'^q^0 , . 

q{x,t = 0) = qo{x) 

where qoix) —^ sufficiently rapidly as \x\ oo. Just as KdV is associated with the 
Schrodinger operator, NLS is associated with a first order, two-by-two scattering 
problem 

d'>P -Z , f iq\ , /rn: 

— t-a^ip + \ ._ n I V : < a; < oo (58) 



dx 2 \-iq 

where "^3 ^ ^ q '^l ^ ^^"^ third Pauli matrix. Let r — r(z) be the reflection 

coefficient for H58|) with q — qo. The map R : q i-^ r is the analog for NLS of the 
OP maps F and Now, for fixed x and t, let m = m{z; x, t) be the solution of the 
normalized RHP (R, Vx.t) where R is oriented from — oo to -f oo and 

, ,( l-|r(z)|2 re*« \ 

~ \ -fe'« 1 ) , ze R. (59) 

Let rni(x, t) be the residue of to at z = oo, 

mi{x,t) , „ / 1 



m(z; x,t) — I -\ h O , „ 

z 

Then 

q{x,t) = -i(TOi(a;,t))i2 (60) 

How does one prove (|()U|I ? At the functorial level, R is really a map from the 
category of differential operators to the category of RHP's, 

L{q) ^ q^ r t-^ v^.t 

where L{q) = icrs-^ -I- ^ ^_ , and so the key question becomes: "How is the 

differential operator encoded into the formalism of RHP's?" 
To answer this question, observe that 

ip = ip{z;x,t) = m{z;x,t)e'^''^ (61) 

solves the RHP 



• il}{z\x^t) analytic on C\R 

/ 1 - |r(z)P r{z) \ 
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where the jump matrix is now independent of x and t. Differentiating with respect 
to a;, we obtain 

l-|r(z)|2 r{z) 
—r(z) 1 

from which it foUows that T = ipx'4'~^ has no jump across M, and hence is entire. 
But as z ^ oo, 



x + 



T — rrixTn 
. era 



2 



A + O 



for some constant matrix A. By Liouville, we must then have T = iz 



. 0-3 



q 
q 



A or 
(62) 

, and hence 



Simple symmetry considerations imply that A is of the form 

we recover the differential equation (|58|1 . Differentiating the V'-RHP with respect 
to t yields similarly an equation of the form 

V't = (63) 

for some explicit matrix B — B{z,q,qx). Cross-differentiating l|62|l and (j^Sl, 
{4'x)t — {'4't)x, then yields the NLS equation (|57|l . It is in this way in general 
that identities and differential relationships are encoded into the RHP. 

To apply the above methodology to OP's, consider the solution X'"' of the RHP 
(M, v) above. Observing that satisfies the same jump relation as X'") across 

M, we conclude as before that T = x("+i)(X("))~i is entire. But 

r==x("+i)(z)(x(")(z))-i 




and again by Liouville we conclude that 



X("+i)(z) = (^z 



1 




(n+1) 



1 




1 




X^"Mx(")(z) (64) 



from which the three-term recurrence relation © now follows by a simple compu- 



tation. Moreover, if we replace the weight uj(x) with LUt{x) 



^2xt 



u;{x) 



/jje2^*w(s)ds 



W^(")(z;t) = X(")(z;t)e(*"+9(*»'' 



git) 



solves the RHP (K, v) with jump matrix v = 



4 log 

1 uj{x) 
1 



2st 



uj{s)ds , 



then 



(65) 



independent of t. Dif- 



ferentiating with respect to t, we obtain as above a differential equation for 
|w^(") = ((z + 5)^3 + x[-'>a, ~ a3Xj"))M^(") . 
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Using r to denote the shift operator, TVF^") = ^"("+1), equation ^ takes the 
form 

Cross- "differentiating" ((HS) and JHII), — rVF^") = T— — -, one is led immediately 

dt dt 

to the Toda flow 

In another direction, if = e^^'^^\ ^i^) — 7ma;^'" + - • Im > 0, then L/'^"-' = 

X(")e5^*^^^°'^ satisfies a jump relation across R with jump matrix ^ — ^ q J 

which is independent of z, and by the above general methodology this leads to a dif- 

ferential equation for t/*^"' with respect to z, — - — = DU^''^\ for some explicit D. 

dz 

d dU^'""' 

Cross- "differentiation" , —Ft/'-"'' =r — ; — , then leads to so-called "string equa- 

dz ^ dz 
tions" for the recurrence coefhcients a„,5„. 

Applying the above methodology to the RHP {S^,v) for OP's on the unit cir- 
cle, we obtain, in particular, simple and direct proofs of Szego recurrence, Geron- 
imus' Theorem on the Schur iterates, and the Pintcr-Nevai formula (see | Sim2| . 
and below). Indeed, let y*^"-* solve the RHP {S^, v) above. Then one observes that 

Y{n) ^ y(n+i) f q ^ j Satisfies the same jump relation as F^") across 5^, 



— n 



" - 1 

and hence (y("))^i is entire. As before, this leads to an equation for y^") and 
which takes the form 

for suitable constants a„,5„,c„. Furthermore (det y("))+ = (det y*^"))_ det v — 

(det y("))_, and so det F^") is entire. But det = det ^F^") ^ ^ ^ ) ) ^ 

as z — > oo, and hence det F^") = 1. Taking determinants of both sides of Ht)7|) . we 
find the relation 

dn = Kcn ■ (68) 
IFtouli the first column of 1)6 7() we obtain the relations 

$„+i = (z -f a„)$„ - 4-i^r,$,*-i (69) 

-kI^I = c„$„ - kLi*:-i • (70) 
Eliminating we obtain the Szego recurrence relation 

$„+i = z$„ - a„$; (71) 

with 0„ replaced by and with Verblunsky coefficient 

Q-n = kIK ■ (72) 
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Letting z ^ oo in (|70(l we find 

Cn = K^an-l (73) 

and hence by 

a„ = a„a„_i. (74) 
Now consider the second cohimn in (|67|) . Setting 

r„ - C($„L^s-") , i„ = C($>s-"-i) (75) 

and using and (dH), we obtain as in and lt7?H) 

zr„+i = (z + a„Q;„_i)r„ - q;„ ( ^^^^ ) t„_i (76) 
Ehminating as we ehminated above, 1)76(1 and H77|l reduce to 

r-n+l ^Tn- antn (78) 

ztn+i = -a„r„ + t„ . (79) 

Defining 

/„ = (80) 
and using H78|l and (|79|l . we obtain the recurrence relation 

, n>0. (81) 

In particular, for z = 0, we see that 

an = /n(0) (82) 

and so 1)81(1 can be written in the form 

= f-ZJ^^^ , n > 0. (83) 

i-/«(o)/„ ' - 



Finally observe that 



^0 _ J51 s-z'^2ii-is _ JSi s-z 



Mz) = - = - . s = e"> (84) 

S — Z 



-0 4^ 



where d/i(6') = Lu{d) — . 

Geronimus' Theorem (see |Sim2| ') states the following: Let 

F{z)= [ i±^dM^) 
■/si s~ z 

be the Caratheodory function for and let /schur = i f ^^"^ associated 

Schur function. Let {fn)n>a solve the recurrence relation 183() with fn\n=a = /schur- 
Then 

/„(0) = a„ , n > 
where {a„}„>o are the Verblunsky coefficients for dfj.. 
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However, a simple computation shows that /schur is precisely /o in 1)84(1 : Hence, 
using the general methodology for RHP's as above, we have proved Geronimus' 
Theorem. Moreover we have the following formula for the Schur iterates: 



/nW = ^ = ^^V^i;?^ , n>0 (85) 



which reduces simply, using (3.2.52) and (2.2.53) |Sim2j . to Golinskii's formula 

([Siin2,, Thm. 32.7). 

Finally we note from |Sim2| . (1.3.79), together with the simple identity J -jz^dii — 

l + z f that 

/schuA-l - _ (i?n-l - ^A„-l) / ^ - 

where A„_i,_B„_i are the Wall polynomials. But from ((85|l . we obtain 



,"/^^(S(£l£i;^£l(£)£l^)d^(0) + 



^„ ( ^-(-)--.:^-(-)--'' ),d;,(g) + $„(,) / ■ 

Comparing with ((86|l we obtain 

^*Jz) = Bn-i - zAn-i (87) 

or equivalently 

$„(z) = zB:_i - (88) 

which is the Pinter-Nevai formula (see |Sim2| 'l relating the OP's to the Wall poly- 
nomials. 

In addition to the formulae and identities obtained above for OP's using the 
RHP's (R, v) and {S^,v), one can, using RHP's closely related to (M, v) and {S^,v), 
derive formulae for Toeplitz and Hankel determinants, or more precisely "relative" 
Toeplitz and Hankel determinants, that are particularly useful for asymptotic anal- 
ysis. The asymptotic analysis of Toeplitz and Hankel determinants, dating back at 
least to the work of Szego in 1915, is of considerable, and continuing, mathematical 
and physical interest, and we refer the reader to |BW| . |E] and the references therein 
for more information and recent results. The "relative" determinant formulae are 
as follows. 

Let uji{x),uj2{x) > be two weights on M and let Dn{ujiuJ2), Dn{u>2) be the 
Hankel determinants associated with the measures LL!i(x)LL!2{x)dx and LV2{x)dx re- 
spectively. (Here we do not require ujiLj2dx and L02dx to be probability measures.) 
Then 

fog = fdtl Mx) ( i \ogu,{x)\ dx (89) 



Dn{(^2) Jo Jr \dt 

where Ut = 1 — t + tuji{x), < t < 1, and Rt is expressed in terms of the solution 
= ((^t"^^^)jj)i<»j<2 of the RHP {R,vtix)) in ^ with 

.*(.)= (J '^*(-);^(-) 

as follows: 
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Similarly, if uJi{6), uj2{0) > are two weights on , with associated Toeplitz 
determinants A„(cJiW2), A„(cj2) respectively, then 

where L0t{9) = 1 — t + tu!i{9), < t < 1, and Rt{d) is expressed in terms of the 
solution = ((y/"+^^)y)i<jj<2 of the RHP {S\vt{e)) in ^ with 



wt(e) = ( g 1 I. ^ = e 



as follows: 

R.iO) = {{Yr'^)jyr%, - {Y}"^%,{Yr\,) ^ (92) 
where /=^ . 

The functions Rt{x), Rt{0) have the interpretation as 1-point functions 

Rt{x) = {n + 1) dn{x,xi,X2, ■ ■ ■ ,Xn) (93) 



Rt{e) = (n + i) / dtxie,ei,...,0n) (94) 

Jei&S^,l<i<n 

for the random particle ensembles (see |Mehp with distributions 

n 

dfi{xo,xi,...,Xn) = (I/Ztsi) Y]_ {xi-Xj)'^Y\_(^t^^)i^j)d^odxi...dxn (95) 

0<j<k<n j=0 

and 

n 

dni0o, 01, ... , 0„) = (l/Zsi) n 1^''' - I' n('^*^2)(0,)d0od0i . . . dOn (96) 

0<i<j<ri j=0 

where Zr, Z51 are normalization constants. 

Note that on we can set W2 = 1, so that An{u!2) = 1 and (|91|l gives us a 
formula, first derived in ID II, purely for A„(a;i). In the non-compact situation on 
R, this clearly cannot be done and we must always work with relative determinants 
as in 

Formulae are due to Deift |D3j . and may be proved by generalizing the 

proof of (|91|l given in 0^ for the case UJ2 = 1- A key ingredient in the proof is the 
notion of an integrable operator. If S is an oriented contour in C, we say that an 
operator K acting on LP(E), 1 < p < 00, is integrable if it has a kernel of the form 

K(z,z ) = — ; , z,z eij (97) 

z — z' 

for some functions fj,gk G L°°{Y^), 1 < j,k < i. Special examples of integrable 
operators appeared in the 1960's in the work of McCoy, Tracy and others, and 
elements of the general theory were discovered by Sakhnovich in the late 60 's, but 
the full general theory of such operators is due to Its, Izegin, Korepin and Slavnov 
|11KS| in 1990. Integrable operators have many useful properties (see e.g. |D1|1. In 
particular, if K is integrable as in (|97|l above, then so is (1 — K)~^ — 1, 

, Fi(z)G,(z') 
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for suitable Fj,Gk, ^ l£ j, k < I. Furthermore, quite remarkably, the functions 
F — {Fi, . . . , Fi)'^, G = (Gi, . . . Gj)-^ can be computed in terms of the solution 
of a canonical, auxiliary RHP. Indeed, define the jump matrix v = I — 2'Kifg^ 

on S, where / = (/i, . . . , //)^, g — {gi, . . . , gi)'^ , and assume for simplicity that 

e 

fj{z)gj{z) — 0, z £ Y.. Then, if m solves the normalized RHP (E, v), we have 

F = m±f and G = (TO±)~^g. (98) 
The proofs of H89|l and (|91|l proceed by expressing the relative determinants -'^"("i"^) 



^A^^ili")'"' terms of Fredholm determinants of integrable operators 



logdet(l-i^) = / — logdet(l-iA' 



tr ( — - — K I dt 

1-tK 



and then using to express (1 — tK)^^K = ((1 — tK)^^ — l)/t in terms of the 
solution of the auxiliary RHP associated to tK . We shall say more about (|89() and 
()91|l in what follows. 



3. Applications of (M, -w) and {S'^,v): asymptotics 

In this section we consider the asymptotics of OP's, denoted (b) in Section 
2. In Section 2, the goal was to show how a variety of identities, equations and 
formulae, mostly classical and well-known, follow from a single, basic methodology 
in RHP's. Here the goal is to describe new results on the asymptotics of OP's that 
follow from the RH method, utilizing in particular the non-linear, non-commutative, 
steepest descent method introduced in |DZ1| in 1993. Although much was known 
(see [Szep about the detailed asymptotic behavior of classical OP's, like Hermite, 
Laguerre, Jacobi polynomials, etc., both on and off the contour of orthogonality, 
little was known about the detailed asymptotics of OP's with respect to general 
weights. The main tool that makes possible the detailed analysis of the asymptotics 
of classical OP's is the existence of integral representations for these polynomials, 
to which the classical method of steepest descent can be applied (see, for example, 
|Sze| . Section 8.71). For general weights, one may view the RHP's (R, w(a:)) and 
{S^tV{9)) as non-commutative analogs of these integral representations, with the 
non-commutative steepest descent method now playing the role of the classical 
steepest descent method. 

We now describe the steepest descent method for RHP's in broad outline: Un- 
fortunately we do not have sufficient space in this article to describe the method in 
detail. In the case of NLS (cf. (|59|l and JHOJ), we write the solution q{x,t) of the 
Cauchy problem for NLS as a functional /, say, of the data re^^ , 

q{x,t) ^ f{re^'). (99) 

^From (|35|) and (|(jO|l we see that 

f{re'')^(^£f^{s;x,t){cj++u;^)^^^^. (100) 
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Using the factorization 



/ 1 -re'» ^ 






I) 











(cf. (j2Z|l), so that 



we obtain 

q{x,t) = 
For r "small" , we have 

q{x,t) 



'fe-^ 



-re-'" 




27^/12 



(101) 



(102) 



(103) 



((/+ "small" )(w+ +uj^) 



r{s)e 



^("^-*^')— + "small" 
27r 



indicating that the classical steepest descent method can be applied as t —^ oo. 
However, when r is no longer "small" , we see from the non-linear dependence of 
q{x, t) on r in l| 103(1 . and from the matrix nature of the problem, that a non- linear, 
non-commutative version of the steepest descent method is required, and this is 
the kind of method that was introduced in |DZ1| . In the classical steepest descent 
method, the integral localizes as t — > oo to a small neighborhood of the stationary 
phase point(s), 9'{zq) = zq — x/2t in the case of NLS, and an explicit asymptotic 
formula for the solution is then obtained by evaluating a Gaussian integral: in the 
fully non-linear case (see |DZ2j |DZ4j ) the RHP {R,Vx,t) locahzes to a RHP in the 
neighborhood of the stationary phase point zq = x/2t, and an asymptotic form for 
the solution 

1 



(104) 



is then obtained by solving this local RHP explicitly (in terms of parabolic cylin- 
der functions, as it turns out). The asymptotic form 1(1 04() was first obtained by 
Zakharov and Manakov |ZM| . by other means. In situations where there is more 
than one stationary phase point, for example for MKdV, where 9 — xz ~V Atz^ 
with stationary phase points ±zo = ±\/—x/12t, the long-time behavior of solu- 
tions of MKdV (see jPZlj ) is a superposition of NLS-like contributions from +zo 
and —zq, as long as these points remain separated, i.e. ^ > c > 0. However, 
in the space-time region where — a;/12t 0, and hence +zo —zq, one is in a 
non-linear "caustic" region which is manifested by the solution taking the form of 

a self-similar oscillation, q{x,t) ~ - — -yj^u{x/ {3ty^^), where u is a solution of the 

Painleve II equation u"{t) ~ tu + 2v?' (see |DZlp . 

Up till this point, the RH asymptotic theory proceeded as a non- linear analog of 
the classical steepest descent method in which all the phenomena that arose could 
be viewed as non-linear counterparts of phenomena that had already arisen in the 
linear, scalar situation. However, with the analysis of the collisionless shock region 
for KdV (see |DVZlj . [DZ3| ). and the analysis of the asymptotic behavior of solu- 
tions of the Painleve II equation, it began to be clear that there were phenomena 
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Figure 2. The contour S 



inherent in the non-hncar steepest descent method that had no analog in the clas- 
sical situation. Most importantly, it became clear that instead of stationary phase 
points, one could have "stationary phase lines" in which case all the points on some 
interval in C contributed equally to the asymptotic behavior of the solution of the 
problem. Moreover, in place of modulated linear oscillations as in H104() . one would 
now have genuinely non-linear oscillations described in terms of Jacobi's sn and cn 
functions, etc. A systematic extension of the steepest descent method to allow for 
such "stationary phase lines" and genuinely non-linear oscillations was presented by 
Deift, Venakides and Zhou J2VZ2^ in the context of their work on the zero disper- 
sion problem. Soon thereafter, using the methods in |DVZ2| together with recent 
developments in the theory of logarithmic potentials with external fields (see |STj , 
and also (DKM ), the authors in DKMVZ2^ derived so-called Plancherel-Rotach 
asymptotics for OP's with measures of the form 

e-^(^)dx, V{x) = 7a;2™ + fe^m-i ^ . . . ^ > (^qS) 

and in |DKMVZT] . for measures of the form 

e-''Qi^)dx, Q(a;)/log|a;| ^ +00 as |a;| c3o, (106) 

where Q{x) is real analytic on R. As described in |DKMV!Z2] one obtains as rt ^ 
00 precise pointwise asymptotics for the OP's Pn{z) for all 2: e C, as well as 
detailed asymptotics for a„, 5„, 7„ and the zeros of pn{z). In the special case 
^-n{x -tx )^^^ Bleher and Its jBI| obtained asymptotics for the associated OP's 
using RH techniques and a mixture of steepest descent/isomonodromy ideas. 

In broad outline the method proceeds as follows. For weights e~^^^^ as above 
one first scales x xn^^^"^ so that e^^*-^^ ^-nVnix) ^ where Vn{x) — jx^"^ + 
_s_^2m-i . . . ^ Next, one considers the so-called equilibrium measure d^eq for 
the logarithmic potential problem associated with OP's (see |ST|V By |DKM| . for 
weights e-^'^riix) qj. ^-nQix) above, dfj,cq is supported on a finite union of disjoint 
intervals U'l^i{ai, hi), J < 00 (in the case e~"^"(^\ J = 1). Next one introduces the 

so-called "5" function, g{z) = / log(z — s)d^cq{s) ^ logz as z — s- 00. Along with 

Jr 

dficq, the logarithmic potential problem also produces a Lagrange multiplier £, and 
we set = eT^'^^A:(")(z)e""»('^)''3e"'^''^ One observes that now solves a 
normalized RHP (M, v) for some explicit jump matrix v. In the key step, the RHP 
for X'"' is now deformed to a RHP on a contour S of the form shown in Figure 
2. By the properties of g{z), or more properly, the properties of d/ieq, it turns out 
that as n — *■ 00, v, the jump matrix for the deformed RHP on S, converges 

v{z) / (107) 

exponentially for all z G S\ u/^^ [ai,bi\. Thus as n — > 00, the RHP reduces to a 
limiting RHP on the union of intervals U'-^i[ai, h]. On each of the intervals {ai,bi), 
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v{z) has the simple form (^^ i q ^ ^^'^ ^^^^^ hmiting RHP can be solved explicitly 

in terms of the function theory on the hyper-elliptic Riemann surface obtained 
by gluing together two copies of C\ uf^^ i'^iih) in the standard way. However, 
the convergence rate in (|107|l is not uniform, becoming slower and slower as z 
approaches the end points ai,bi. The natural topology for RHP's is convergence for 
the coefficients of v in L^'nL°°(S) (cf. l|3Q(IH33(l 'l. and the lack of uniform convergence 
in Hl()7(l constitutes the major technical difficulty in implementing the steepest 
descent method as described above. We refer the reader to IDKMVZTl IDKMVZ'2] 
for more details. 

We now consider the relative determinant formulae H89(l and (|91|) and their as- 
sociated RHP's (R,Dt(x)) and {S-^,Vt{0)) respectively. The celebrated strong Szego 
limit theorem, in the definitive form due to Ibragimov (see jSim2j for many proofs 

do 

and much historical discussion) states that if dfj,{6) — s^^'^^'^t;—, and V{9) has 

ZTT 

oo 

Fourier coefficients {14} satisfying A:|VfcP < oo, then as n — > oo 

fc=i 

C30 

lnA„(e-^) = (7i + l)fo + ^fc|14p + o(l) . (108) 

fc=i 

In addition to the many proofs in |Sim2| . (jl()8|l can also be proved, under cer- 
tain additional smoothness assumptions on V{9), by applying the steepest descent 
method to the RHP's {S^,vt{0)), < t < 1. The situation is simpler than in 
[DKMVZI. DKMV!2S] , but the argument in this situation is particularly illustra- 
tive of the emergence of a "stationary phase line" : details are given in jni]. There 
is also a version of the strong Szego limit theorem for block Tocplitz determinants 
(see |W1| . |W2| . and also (Bot, for more recent results). In the block Toeplitz case, 
the analog of (|108|l contains a certain Fredholm determinant which is difficult to 
evaluate in elementary terms. In certain cases the method in Dl extends to the 
block Toeplitz case and, quite surprisingly, the term corresponding to this Fredholm 
determinant is evaluated automatically (see |LTK| 'l. 

In |BCW1| (see also |BCW2| ') the authors state the following analog of the Szego 

2 

strong limit theorem for the case of Hankel matrices. Let W2 = and let 

u)i{x) > have the property that lui{x) — > 1 sufficiently rapidly as — > oo. Then 
as n ^ oo 



ln%^^^^^ / logc.,(.)dx + -L / |fc| \m\^dk + oil) (109) 

where f{k) = j_ / {logLUi{x))e^'^''^ dx. Using (|S^ . this result can also be proved 
V 27r J 



( |D4j ) using the steepest descent method, not only for lu2 — e^^ , but also for more 
general weights, uj2 = e"^*^^\ ^ix) = 72;^™ -I- . . ., 7 > 0, as above. 

Riemann-Hilbert techniques and the RH method are useful not only for asymp- 
totic evaluation, but also for estimation. For example, let o-!{0) E L°^{S^) be a 

e^^^^uj(6) — , keZ. Let 
27r 

{{T{uj))jk)'^k=o — i^j-k)'j°k=o denote the Toeplitz matrix associated with uj acting 



RIEMANN-HILBERT METHODS IN THE THEORY OF ORTHOGONAL POLYNOMIALS 19 



on£+ = lu = (mo,mi,. •■) : l^^l^ <°o\^ and let ((T„(w))j7c)"fe=o = K-fe)"fe=o 
denote the leading {n + 1) x (n + 1) section of T{uj). If w is in the Wiener space 



Icjjl < oo > with uj{9) > 0, then, by a well known theorem of 

j=-oo 

Krein, {T{llj))^^ exists as a bounded operator in The question is the following: 
How closely does {Tn{uj))~^ approximate {T{uj))~^ for n large? Let i' = {i^k)kez be 
a Beurling weight (see e.g. jSim2 |): Thus i^k ^ ^, I'k = ^-k and Vj-\-k < VjVk for all 
j, fc e Z. In particular, ((1 + |fc|)^)fcez, £ > 0, and (e"l'°l)fegz, a > 0, are Beurhng 
weights. Define the Beurling class 



W,^LeL\S'):Y.v,\u:,\<^\. 



Clearly C W° for any Beurling weight v. Let G 14^,^ for some v and assume in 
addition, for simplicity, that the weights increase on i.e. Vj < Vk for < j < fc. 
Then the following is true |DO| : for n sufficiently large and < j, k < n, 

{Tn{uj))]k - {n^))]k\ < c.(^) min(i.-^,_„ (110) 

for some constant c,y{uj). Thus for < j,k < n, {Tn{ijj)') is a good approximation 
to (r(a'))^.J^, apart from the lower right corner j ^ k ^ n. This estimate is a 



generalization of an earlier estimate due essentially to Widom (see jBS| for refer- 
ences and further discussion). The proof of (|ll()|l in jPOj uses RH techniques in 
an essential way closely related to the proof of JHSJ. The paper also contains other 
results for orthogonal polynomials on the unit circle, including a new RH proof of 
the reverse statement in Baxter's theorem (cf. |Sim2| '). Interestingly, the Borodin- 
Okounkov operator jBO| . or more properly, the Borodin-Okounkov-Case-Geronimo 
operator, which has emerged recently as a powerful tool in the analysis of Toeplitz 
determinants, arises naturally in the analysis in |DO| . 

The steepest descent method for varying weights uj{x) = e~"'3(a:) jjj |DKMV!ZT] 
can also be applied to orthogonal polynomials on the unit circle with varying weights 
uj{0) — e^"'^'^^). For example in their analysis of the length In — /^(ti') of the longest 
increasing subsequence of a random permutation tt on n letters, the authors in jBD,T| 
prove that 

'L, - 2x/n 



lim Prob " ^ J <t] ^F2{t) (111) 

where F2{t) is the Tracy- Widom distribution function for the largest eigenvalue 
of a random matrix from the Gaussian Unitary Ensemble. The proof of (|lll|l in 
|BD.Tj reduces, by a formula of Gessel, to the analysis of the Toeplitz determinant 

A„_i(e'"^°'^^) where s = (n + 1) ( 1 , — I as n — > cxd, and where t is 

\ 2^/'^(n -t- 1)"^'""/ 

the same as in (|lllf) . As indicated above, the method of BDJ is modeled on the 

RH steepest descent method in |DKMVZlj . The same RH problem with weight 

gscose ^1 ^jgQ appears in the work of Baik and Rains (BR| in their analysis of 

monotone subsequences of involutions. 

The steepest descent method for OP's {0n} on the unit circle can also be used 

to obtain detailed information on the zeroes of the (^„'s as n oo (see |MMSj ). 
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In a further development, the authors of jMMj have introduced an extension of 
the steepest descent method to non-analytic weights, obtaining in particular new 
results for the zeros of OP's on the unit circle for such weights. 

Throughout this paper we have restricted our attention to measures that are 
smooth as in (|37|l and H47() . The OP problem for general measures is then ana- 
lyzed (cf. Technical Remark above) by approximating the measure appropriately 
by smooth measures d^f , and then taking the limit as e — s- 0. This approach works 
well for the derivation of equations, formulae, etc., but for asymptotic questions 
one clearly needs a different approach. Recently remarkable connections have been 
discovered (fj') between various combinatorial problems - random growth mod- 
els, random word problems, tiling problems - and certain polynomials orthogonal 
with respect to discrete measures. The polynomials that arise include the classical 
Meixner, Charlier, Krawtchouk and Hahn polynomials (sec Sze ). Related discover- 
ies have also been made in the representation of the infinite dimensional symmetric 
and unitary groups |B01|[5n2) . The Meixner, Charlier and Krawtchouk polyno- 
mials all have convenient integral representations (see Pze] *) and their asymptotic 
behavior can be read off using the classical method of steepest descent. This is 
unfortunately not the case for the Hahn polynomials (such polynomials are needed 
in particular to describe the tiling of hexagons by rhombi). It turns out, however, 
that discrete OP problems can be rephrased in terms of a discrete RHP, which is an 
analogue of the continuous case, and which was introduced by Borodin, along with 
a theory of discrete integrable operators, in 0. In a significant further develop- 
ment of the nonlinear steepest descent method, the authors in jBKMMj extended 
the method to a wide class of discrete RHP's which includes the discrete RHP for 
the Hahn polynomials (as well as the other three discrete OP systems mentioned 
above). The relevant limit here is when the order of the OP's p„ becomes large 
and simultaneously the spacing between the points in the measures goes to zero 
at a prescribed rate (see jBKMMj ). In this way the authors are able to analyze 
the Hahn polynomials asymptotically, proving en route a conjecture of Johans- 
son in U' that for hexagonal tiling the so-called "arctic circle" of jULPj exhibits 
Tracy- Widom fluctuations as in (jllljl above. In jB02j the authors also consider an 
asymptotic problem for Hahn polynomials using a discrete RHP, but the relevant 
limit is different from that in jBKMMj . 

Many researchers are currently involved in the application of RH techniques 
to the theory of OP's. In addition to those mentioned above, the list includes 
Chen, Claeys, Kapaev, Kitaev, Kuijlaars, van Assche and Vanlessen, amongst many 
others. Because of space limitations, however, we unfortunately cannot describe 
their work in any detail, and we must refer the reader to the literature. 

4. Related areas 

In this final section we will describe, very briefly, various areas related to OP's 
in which the RH method plays a role. 

We first consider random matrix theory (RMT), which has been a major source 
of questions and challenges to OP theorists for over 40 years (see e.g. jMehj and 
|U2j '). The situation is as follows. A Unitary Ensemble (UE) is an ensemble of 
N X N Hermitian matrices {M = M*} with probability distribution 

PN{M)dM = ^e-*'^(*-^)dM (112) 
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where 

• dM denotes Lebesgue measure on the algebraicahy independent elements 
of M. 

• W{x) is a real- valued function that goes to +00 as |a;| 00. The case 
W{x) = gives rise to the Gaussian Unitary Ensemble (GUE). 

• Zn is a normalization coefficient. 

"Unitary" refers to the fact that the distribution H112|l is invariant under unitary 
conjugation, M — > UMU*, U unitary. The Universality Conjecture for UE's (see 
|Mehj and p32 ) states, in particular, the following: Given W, if Jm — C]y+S]y{—t, t) 
is a suitably centered and scaled interval in R, then as ^ 00, P{Jn) — Prob(M : 
Mhas no eigenvalues in Jjv) converges to a universal limit independent of W, 

lim P{Jn) = det (1 - St) (113) 

where St is the trace class operator with kernel St{x,y) — ^^'^j!^^y^'^ acting in 
L^{—t, t). The specific form of the weight e~^'^^^ dx is reflected only in the precise 
values of Cat and sat. OP's enter the picture because of the celebrated result of 
Gaudin and Mehta (see |Mehj ) that if i? C M is a Borel set, then 

Prob(M : Mhas no eigenvalues inB) = det(l - Kn,b) (114) 
where ifjv.B is the finite rank operator with kernel 

N~l 

/^^(a:,y)= ^p,(x)p,(y)e-^'^(-)e-^^(^) (115) 
j=o 

acting on L^{B), and {pj}j>o a-re the orthonormal polynomials Q with respect 
to the weight e"^^^-' dx. Hence the question of proving universality as in H113|) 
becomes a question of deriving the appropriate asymptotics for OP's, and this is 
the main scientific content of |UKMVZ2j . |UKMVZ1| . [DS and |BI|. Of course, 
if the weight e~^'^'^^ dx is classical, e.g. W{x) = x^, and the asymptotics of the 
associated polynomials {Pj}j>o can be derived from an integral representation, 
then universality for these ensembles can be proved without recourse to the RH 
steepest descent method, and this has been done by various authors (see jI)KMVZ2], 
[HKMVZI for references to the literature). 

Orthogonal ensembles (OE's) oi N x N real symmetric matrices {M = M = 
M^} and Symplectic Ensembles (SE's) of 2N x 2N Hermitian self-dual matrices 

{Af — M*,JMJ'^ — Af'^}, where J — diag(r, . . . , t), t — , equipped 

with invariant weights analogous to H112() . are more difficult to analyze. Firstly, in 
the place of determinantal expressions as in H114(l , one obtains Pfaffians (see |Meh| 
for classical ensembles, jTWj for the general case) 

Prob(M : Mhas no eigenvalues inB) = (det(l - Kn^b))^^^, (116) 

and, moreover, the operators i^Ar.B are now 2x2 matrix operators with kernels 
(^K]sf^ij{x, y)) j^<;j ^<2> x,y G B. In contrast to Ijll5|) . these kernels are most naturally 
expressed in terms of certain skew-orthogonal polynomials (see jMehj ^. but for gen- 
eral weights e"'^^^-' dx the asymptotic behavior of such polynomials is not known. 
However Widom jW8j has shown that if W /W is rational, then (^KN^ij{x,y)) can 
be expressed conveniently in terms of the orthonormal polynomials {pj}j>o with 
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respect to the weight e^^'^^^ dx, so again, as in the unitary case, the question of 
universaUty of OE's and SE's becomes a question of analyzing the asymptotic be- 
havior of OP's. The expressions for {KM,ij{x, y)) are now more cumbersome 
than H115|l and significant new technical issues arise, but nevertheless, using the 
asymptotic analysis in DKMVZ2j as a basic ingredient, it is indeed possible to use 
Widom's formulae in |W3j to prove universality for OE's and SE's with weights of 
the form e'^^^^ dx, V{x) = jx'^"' + • • • , 7 > 0. This is the content of [PGl] and 

Biorthogonal polynomials TTk{x) — + ■ ■ ■ , a j{y) — ■ ■ ■ , k,j>0, 

nkix)c7j{y)e-^^^^~^^y^+^^''ydxdy^0 iij^k, (117) 

arise in the analysis of the theory of coupled random matrices. Here V{x) and 
W{y) grow sufficiently rapidly as |a;|, |y| 00, and t ^ 0. Various RH problems 
have been proposed to analyze these polynomials (see, in particular, |BEH', 
|KM| and the references therein), but the analysis of the asymptotic behavior of 
these RHP's is still at a preliminary stage. 

For TO > 2, let n = (rii, n2, . . . , Um) be a vector of non-negative integers, and 
let iOi{x) > 0, . . . ,u!mix) > be weights on M with finite moments. Let \n\ = 
ni + - ■ ■ + nm- Multiple orthogonal polynomials (see |Apt|) of type I are polynomials 



(k) (k) 

An for k = 1,2, ... ,m, deg A„ < rifc — 1 such that the function 



H„ix)^y^Ai^\x)LOk{x) 

satisfies 



k=l 



Multiple orthogonal polynomials Ln{x) of type II are monic polynomials of degree 
\n\ satisfying 

Ln{x)x^u}j{x) dx = for fc = 0, . . . , — 1, j = 1, . . . , to. (119) 

Multiple orthogonal polynomials were first introduced by Hermite in his proof of the 
transcendence of e. In 2000, van Assche, Geronimo and Kuijlaars |vAGK| showed 
that multiple orthogonal polynomial problems of types I and II could be rephrased 
as RHP's analogous to the RHP of Fokas, Its and Kitaev for ordinary OP's, and 
they used these RHP's to derive various properties and relations for the multiple 
OP's. In the last year or two significant progress has been made in extending and 
applying the steepest descent method to RHP's which arise from multiple OP's 
in special cases. We mention, in particular, |BK| . jABKj and |K VaWj . |KSVaW] 
and the references therein: In the first two papers the authors consider a random 
matrix ensemble PN{M)dM = ^^e^^^'^s^-^'-^^^) dAf , with external source A, 
first analyzed by Pastur, Brezin-Hikami, and later by Zinn- Justin. Under certain 
conditions on A, they show that the ensemble can be analyzed as iV — > 00 in 
terms of a 3 x 3 RHP to which an extension of the nonlinear steepest descent 
method can be applied: A new phenomenon now occurs in the analysis, which 
the authors term a "global opening of lenses" (see |ABK| ). In the second two 
papers the authors analyze type I and type II Hermite-Pade approximations to the 
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exponential function, which they are again able to control by applying an extension 
of the steepest descent to a 3 x 3 RHP. 

Riemann-Hilbert techniques can also be used to analyze the asymptotics of so- 
called orthogonal Laurent polynomials. Such polynomials arise in the following way. 
Let V{x) be a real-analytic function on M \ {0} with the property 

V V{x) V{x) 
Ix|^ooln|2:| |a:|^o ln(|a::| '-) 

Orthogonalization of the ordered basis {1, z~^, z, 2:^, . . .} with respect to the 
pairing (/, 5) 1-^ J^f{s)g{s)e~^'^^'^^ds leads to the even degree and odd degree 

orthonormal Laurent polynomials {(f>m}m>o'- 4>2n{z) = ^l^^'^'z^" + • • • -I- 

ei'"^ > 0, - ^ + e:^^ > 0. Recently, 

McLaughlin, Vartanian and Zhou (see |MVZj and the references therein) have used 
RHP-steepest descent methods to analyze the asymptotic behavior of the Laurent 
polynomials <t)2niz),<t)2n+i{z) and their associated norming constants ^n"\S,'"^n-i'^ 
in the limit as 00, N/n 1. The work of McLaughlin et al. involves 

significant extensions of the steepest descent method: Such extensions are needed in 
order to overcome the new difficulties introduced into the problem by the singularity 
of the potential V{x) at x = 0. 

Finally, there are problems in which the asymptotic behavior of the system at 
hand is described by OP's. This happens, in particular, in the case of the so-called 
Toda rarefaction problem (see [DKKZp . Here one considers the initial-boundary 
value problem for the Toda lattice 

= e^"-^-^" - e^"-^"+\ n>l (120) 

where for some a > 

Xn{0) = an, n>l; 

i„(0) = 0, n>l, ^'^'> 

and the driving particle moves with a fixed velocity 2a 

xo{t) = 2at, t > 0. (122) 

Making the change of variables x„ — > an + y„ one sees that, apart from rescaling 
time, one can always assume without loss of generality that a = in H121|l . One 
may think of H12U|) - (|122|I as a cylinder of particles {xn}n>i driven by a piston xq- 
If a > 0, one has the (Toda) shock problem f |VDOp and if a < one has the 
(Toda) rarefaction problem. In the rarefaction problem, if \a\ is sufficiently large 
(|a| > 1 turns out to be the critical region) one expects that the piston will separate 
from the "gas" {xn}n>i and cavitation will occur. This is indeed what happens: 
if a < —1, the authors in [DKKZj show, using the RH steepest descent method, 
that as t — *■ 00, the solution of the Toda lattice splits into two parts, I-I-II. Part I 
models the cavitation and Part II is an exponentially decreasing error term. Quite 
remarkably. Part I is constructed from the solution of an associated OP problem, 
which turns out to be the Fokas, Its, Kitaev RHP in disguise. We refer the reader 
to ,DKKZ| for details. 
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